Conduction Mechanism in a Molecular Hydrogen Contact 
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We present first principles calculations for the conductance of a hydrogen molecule bridging a pair 
of Pt electrodes. The transmission function has a wide plateau with T fa 1 which extends across the 
Fermi level and indicates the existence of a single, robust conductance channel with nearly perfect 
transmission. Through a detailed Wannier function analysis we show that the H2 bonding state is 
not involved in the transport and that the plateau forms due to strong hybridization between the 
H2 anti-bonding state and states on the adjacent Pt atoms. The Wannier functions furthermore 
allow us to derive a resonant-level model for the system with all parameters determined from the 
fully self-consistent Kohn-Sham Hamiltonian. 

PACS numbers: 73.63.Rt,73.20.Hb,73.40.Gk 



The study of electron transport through single 
molecules has evolved during the last decade as new ex- 
perimental techniques have made it possible to produce 
atomic-scale contacts with a few or even a single molecule 
suspended between macroscopic electrodes 0, 0, 0> 0- 
At the same time theoretical efforts have been made 
to describe and understand the experiments from first 
principles 0, 0, • The connection between experiment 
and theory, however, has been complicated by the cru- 
cial but in practice uncontrollable atomistic details of 
the contact region between the molecule and the leads. 
While the majority of previously investigated molecules 
have shown a conductance much lower than the quan- 
tum unit, Go = 2e /h, Smit. et al. recently measured a 
conductance close to IGq for a hydrogen molecule bridg- 
ing a pair of Pt electrodes |S|- The result immediately 
raises the question: how can a hydrogen molecule which 
has a closed shell configuration and a large energy gap 
be conducting? Despite the simplicity of the system, 
there are still considerable disagreements among the re- 
ported calculations for the conductance of the hydro- 
gen bridge. Quantitatively, values of 0.9Gq H, U and 



(0.2 — 0.5)Go 10] have been published by different group 
using similar methods. Perhaps even more importantly, 
the physical explanations for the obtained results are very 
different. Indeed, both the bonding 0,0] as well as the 
anti-bonding Q state of the H2 molecule have been pro- 
posed as the current-carrying state. 

In this letter we present conductance calculations 
based on Density Functional Theory (DFT) showing that 
a hydrogen molecule bridging a pair of Pt contacts can 
have a conductance close to IGo and we explain the 
physical mechanism behind this result. The transmis- 
sion function is found to have a characteristic plateau 
with T ps 1 in an energy window of 4 eV around the 
Fermi level, indicating the existence of a single, very ro- 
bust conductance channel with nearly perfect transmis- 
sion. By performing a Wannier function (WF) analysis 



we can directly study the transmission through the H2 
bonding and anti-bonding states separately. The results 
clearly demonstrate that the bonding state takes almost 
no part in the transport and that the plateau is a result 
of a strong hybridization between the H2 anti-bonding 
state and a combination of d- and s-like orbitals located 
on the neighboring Pt atoms. The analysis furthermore 
allows us to determine characteristic model parameters 
from first principles which in turn provides a very simple 
description of the system. 

To describe the molecular contact we use the supercell 
shown in the inset of Fig. It contains the H2 molecule 
anchored between two 4-atom Pt pyramids which again 
are attached to Pt(lll) surfaces We calculate 

the conductance of the relaxed structures assuming that 
the electrons move phase coherently through the con- 
tact and are influenced only by the self-consistent Kohn- 
Sham potential. In this case the conductance is given 
by G = GqT(sf), where T(sf) is the transmission func- 
tion at the Fermi level jl5| . The transmission function 
is found usin g th e Green's function method described in 
Refs. [la, Il7l Il8j . In this approach the system is divided 
into three regions: A left lead, L, a right lead, R, and a 
central region, C. The leads are assumed to be periodic 
such that all scattering takes place in G. In our case G 
coincides with the supercell of the DFT calculation and 
the leads are bulk Pt(fcc) described in a supercell con- 
taining 3x3 atoms in the transverse plane to match the 
central region at the interfaces. The transmission func- 
tion is then given by the formula 



T{e) = Tr[G r c (e)T L (e)G a c (e)r R (e)], 
j is 

tcring region 



(1) 



where G r c is the retarded Green's function of the scat- 



G r c (e) = [(e + ir,+)S - £ L (e) - E R (e) - H c ] 



(2) 



Here He and S are the Hamiltonian and overlap matrices 
of the central region, rj + is a positive infinitesimal and S Q 



2 



is the self-energy from lead a. The coupling strength of 
lead a is given by T a = i(S a — Sjj. 

We use partly occupied Wannier functions, {^ no }, (see 
below) as basis functions in each of the three regions 
(a = L,R,C). Due to the limited size of the super- 
cell in the plane perpendicular to the transport direction 
the conductance should be calculated as an integral over 
the Brillouin zone in the corresponding plane. We thus 
form the Bloch states ipk ±na ( r ) = ^ R e lU± ' R± (f> na (r — 
R_l), where runs over supercells in the transverse 
plane. For each we obtain a Hamiltonian matrix 
H(k±) n a,mp = (k±na\H\k±mf3) which in turn leads to 
a conductance G(k±) through Eq. QJ. The integrated 
conductance can then be approximated by the finite sum 
J2 k± w(kj_)G(kj_), where w(k±) are appropriate weight 
factors. 

We focus on a single, fully relaxed contact character- 
ized by the bond lengths dn-H = 1.0 A and dp t _H = 
1.76 A. The vibrational modes of the hydrogen molecule 
in this configuration are in fair agreement with new ex- 
perimental results 20] . In Fig. |I]we show the trans- 
mission function calculated using 8 irreducible k-points 
to sample the transverse BZ. The same curve calculated 
within the widely used r-point approximation is shown 
for comparison. The two curves have essentially the same 
features, however, the r-point curve has more structure. 
This is because the k-point sampling provides the correct 
smearing of the electronic structure in the leads which ef- 
fectively washes out features related to single points in 
the transverse plane of the lead Brillouin zone. An in- 



teresting feature of the transmission function is the wide 
plateau with T«l extending across the Fermi level. We 
shall refer to this plateau as the lGo-plateau. 

To gain insight into the formation of the lGo-plateau 
we perform a Wannier function analysis. The WFs are 
defined as linear combinations of the Kohn-Sham eigen- 
states with the expansion coefficients chosen to make the 
WFs orthogonal and maximally localized. By including 
selected unoccupied eigenstates in this construction we 
can obtain good localization properties of the WFs also 
for metallic systems [^, |2^- We stress that the mini- 
mal WF basis set retains the accuracy of the plane wave 
DFT calculation since the WFs by construction span the 
eigenstates below a certain energy which has been set 
to 4 eV above the Fermi level in the present calcula- 
tion. The transformation results in the following set of 
WFs: For each Pt we obtain 5 <i-orbitals centered at 
the atom and a single cr-orbital located at an intersti- 
tial site. For each hydrogen we find an s-orbital, \i) 
(i = 1,2), which is slightly elongated towards the con- 
tacting Pt atom. We proceed by transforming the hy- 
drogen s-orbitals into bonding and anti-bonding combi- 
nations \b) = (|1) + |2))/a/2 and \a) = (|1) - |2))/v / 2. 
| a) and \b) are the only states with significant weight 
on the molecule and provide two conductance chan- 
nels well separated in energy. The on-site energies are 
(b\H\b) = -6.4 eV and (a\H\a) = 0.1 eV relative to 
the Fermi level of the metal. By cutting all coupling 
matrix elements involving the bonding, respectively, the 
anti-bonding state we can directly test their individual 
contributions to the total conductance when interference 
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FIG. 1: Calculated transmission for the molecular hydrogen 
contact shown in the inset. For comparison both the k-point 
sampled transmission and the F-point transmission are shown. 
The r-point transmission has more structure, however, the 
qualitative features of the curves are essentially the same. 
The wide plateau with T m 1 extending across the Fermi 
level indicates a single, very robust conductance channel with 
nearly perfect transmission. 
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FIG. 2: Calculated transmission for the structure shown in 
the inset of Fig. when all coupling to the bonding, respec- 
tively, the anti-bonding Eb state has been cut. The full trans- 
mission has been repeated for comparison. The narrow peak 
around — 7 eV is clearly due to the bonding state, while the 
peak at —2 eV and the wide plateau around the Fermi level 
almost exclusively involve the anti-bonding state. 
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is neglected. The result is shown in Fig. [21 The nar- 
row peak just below —7 eV is completely gone when the 
bonding state is removed but is not affected by the ab- 
sence of the anti-bonding state. The peak is thus clearly 
due to transmission through the bonding channel which 
is in good agreement with the calculated on-site energy 
of | b). In the energy regime -6 - -4 eV both the bonding 
and anti-bonding states contribute to the transmission. 
For energies above -3 eV the removal of the bonding state 
has little effect on the transmission which must therefore 
be ascribed to the anti-bonding state. A small exception 
to this is the narrow peak at -1 eV which is caused by hy- 
bridization of the bonding state with Pt d z 2 -orbitals on 
the contacting atoms. Overall we can conclude that the 
peak at -2 eV and the lGn-plateau which determines the 
conductance are due to transmission through the anti- 
bonding state. 

The fact that the bonding state takes almost no part in 
the transmission around the Fermi level allows us to de- 
scribe the contact by a resonant level model |2^| with all 
parameters determined from first principles. In the reso- 
nant level model we consider a single level, \a), of energy 
£a = ( a \H | a) coupled to infinite leads via the matrix ele- 
ments tk a = (ka\H\a), where {|&a)} is a basis of lead a. 
The model has served as starting point for many more ad- 
vanced studies such as shot noise, electron-electron and 
electr on-p honon interactions in resonant tunneling sys- 
tems [3 EE HI3 ■ A particularly useful formulation of 
the model can be obtained if we introduce the group or- 
bital of lead a by \g a ) = c a P a H\a) , where P a is the 
orthogonal projection onto lead a and c Q is a normal- 
ization constant. By supplementing the group orbital by 
orthonormal states {|fca)} we obtain a new basis with the 
key property (ka\H\a) = for all k. The level is thus 
coupled to the lead via the group orbital only. Since the 



rj H 2 anti-bonding state \) 



%9* <*£ 

' J Group orbitals \ 



FIG. 3: Contour plots of the orbitals determining the trans- 
port properties of the hydrogen contact: the H2 anti-bonding 
state, \a), and the corresponding left and right group orbitals, 
\gi), \9r)- The left (right) group orbital has been constructed 
by applying the DFT Hamiltonian to \a) and then project- 
ing onto the Wannier functions of the left (right) part of the 
contact. 



contact is symmetric (gi\H\a) — (gn\H\a) = V. The 
imaginary part of the level self-energy, A a , is directly re- 
lated to the DOS of the group orbitals calculated with 
V = 0: A a = i:\V\ 2 {pl + p° R ) = {T L + T R )/4, The real 
part of the self-energy is the Hilbert transform of A Q . 
For a symmetric contact we have p L = p R = p a g and the 
transmission in Eq. Q takes the simple form 



T{e) = 2^\V\ 2 Pq (e) Pa (e) 



(3) 



Since p a can be obtained from S a and e a this expression 
shows that the transmission is determined by the three 
quantities p g , V and e a . 

By applying the DFT Hamiltonian to \a) in the WF ba- 
sis we construct the group orbitals of the H2 anti-bonding 
state. Contour plots of the orbitals are shown in Fig. [3] 
The group orbital is mainly composed of the d z 2 -orbital 
of the apex Pt atom and three cr-orbitals centered within 
the Pt pyramid. We calculate p° g for the uncoupled sys- 
tem by cutting all coupling matrix elements to \a). The 
result is shown in the upper panel of Fig. ^together with 
p a and the transmission function. The pronounced peak 
at -1 eV is due to the d z 2 -orbitals on the apex Pt atoms. 

If we neglect the narrow peak at —1 eV, p° g can be 
described by a semi-elliptical band on top of a flat back- 
ground, see lower panel of Fig. 0] The coupling and level 
energy can be directly read off the Hamiltonian matrix 
and we find V = 1.9 eV and e a — 0.1 eV relative to the 
Fermi level. It should be noticed that the coupling which 
is relevant for the adsorption of the hydrogen molecule 
to the contact is \[2V = 2.7 eV since the level is coupled 
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FIG. 4: Upper panel shows the transmission together with 
the projected density of states (DOS) for the H2 anti-bonding 
state and the corresponding left group orbital. The lower 
panel shows the same quantities obtained from the single-level 
model when p° g is approximated by a semi-elliptical band and 
we use the coupling, V, and on-site energy, e a , from the first 
principles calculation. 
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by V to both leads. From these parameters we can de- 
termine the Green's function for the level which in turn 
yields p a and T. The result is summarized in the lower 
panel of Fig. 0] Based on the good agreement with the 
first principles results we conclude that the simple model 
indeed gives a realistic description of the system. It is 
then clear that the peak at — 2 eV represents the bonding 
combination between \a) and the Pt band and that the 
IGo-plateau forms because: (i) e a lies close to the Fermi 
level and well inside the relevant Pt band as defined by 
the group orbital, (ii) the width of the renormalized level 
(A a ) is comparable to the band width, i.e. the limit of 
strong chemisorption [23j |. 

The crucial point in the proposed mechanism is the 
strong hybridization of the H2 anti-bonding state with 
the Pt bands around the Fermi level. This picture 
agrees well with the conventional understanding of hydro- 
gen dissociation on simple and transition metal surfaces 
which has been established on the basis of DFT calcula- 
tions [23, ■ The bonding and antibonding states of a 
hydrogen molecule at a simple metal surface are broad- 
ened and furthermore shifted down due to the hybridiza- 
tion with the metal s- and p-states. During the dissocia- 
tion process the antibonding resonance crosses the Fermi 
level and becomes gradually filled with the result that 
the hydrogen-hydrogen bond is weakened. For the tran- 
sition metal the same general picture applies, but the hy- 
bridization with the d-states further affects the antibond- 
ing resonance. The fact that the antibonding state in the 
calculations for the bridging hydrogen molecule between 
Pt contacts is close to the Fermi level is thus an indica- 
tion that the hydrogen- hydrogen bond is weakened by the 
coupling to the metal in agreement with the resulting in- 
crease of the hydrogen-hydrogen bond length. The values 
we find for the positions of the bonding and antibonding 
molecular levels, = —6.4 eV, e a = 0.1 eV are in fact 
quite close to the ones used by Hammer and N0rskov 
(st, = —7 eV, e a = 1 eV) to describe hydrogen in the 
dissociative transition state on metals. This is in clear 
contrast to the studies by Cuevas et al. || and Garcia et 
al. [lfl ] who consider the hydrogen molecule in the bridg- 
ing position to have almost the same bond legnth as 
the free molecule and who report very large bonding- 
antibonding splittings of 23-24 eV which even exceed the 
DFT-PW91 value of 10.4 eV for a free molecule. 

In summary, we have presented first principles conduc- 
tance calculations showing that a hydrogen molecule sus- 
pended between Pt contacts can have a conductance close 
to lGo- Through a detailed Wannier function analysis we 
have identified the conduction mechanism as being due 
to a strong hybridization between the H2 anti-bonding 
state and certain Pt bands. A resonant level model with 
all parameters determined from the self-consistent DFT 
Hamiltonian was shown to account for the important fea- 
tures of the first principles transmission function. 

We would like to thank J. van Ruitenbeek and D. Dju- 



kic for many inspiring discussions. We acknowledge sup- 
port from the Danish Center for Scientific Computing 
through Grant No. HDW-1101-05. 



[1] C. Joachim, J. K. Gimzewski, R. R. Schlittler and 

C. Chavy Phys. Rev. Lett. 74, 2102 (1995). 
[2] M. A. Reed, C. Zhou, C. J. Muller, T. P. Burgin and 

J. M. Tour Science 278, 252 (1997). 
[3] J. Reichert, R. Ochs, D. Beckman, H. B. Weber, M. 

Mayor and H. v. Lohneysen Phys. Rev. Lett. 88, 176804 

(2002). 

[4] B. Xu and N. J. Tao Science 301, 1221 (2003). 

[5] J. Taylor, M. Brandbyge and K. Stokbro Phys. Rev. Lett. 

89, 138301 (2002). 
[6] P. S. Damle, A. W. Ghosh and S. Datta Phys. Rev. B 

64, 201403 (2001). 

[7] M. Di Ventra, S. T. Pantelides and N. D. Lang Phys. 

Rev. Lett. 88, 046801 (2002). 
[8] R. H. M. Smit, Y. Noat, C. Untiedt, N. D. Lang, 
M. C. van Hemert and J. M. van Ruitenbeek, Nature 
419, 906 (2002). 
[9] J. C. Cuevas, J. Heurich, F. Pauly, W. Wenzel, and G. 
Schon, Nanotechnology, 14, R29 (2003). 

[10] Y. Garcia, J. J. Palacios, E. SanFabian, J. A. Verges, 
A. J. Perez- Jimenez and E. Louis, Phys. Rev. B 69, 
041402(R) (2004). 

[11] The surface is modeled as a slab with a thickness of 4 
atomic layers and the supercell contains 3x3 atoms in 
the surface plane. Both the Pt pyramids and the H atoms 
have been relaxed using a plane wave DFT code [l2^ to 
obtain the most stable contact geometry. We use an en- 
ergy cut-off of 25 Ry for the plane wave expansion and 
describe the ion cores by ultrasoft pseudopotentials [l3|. 
To treat exchange and correlation we use the PW91 func- 
tional [l4|. The Brillouin zone is sampled by a single k- 
point along the contact axis (z-axis) and 4x4 k-points 
in the transverse plane. 

[12] B. Hammer, L.B. Hansen, and J.K. N0rskov, Phys. Rev. 
B 59, 7413 (1999); SR. Bahn and K.W. Jacobsen, Comp. 
Sci. Eng. 4, 56 (2002); The Dacapo code can be down- 
loaded at http://www.fysik.dtu.dk/campos 

[13] D. Vanderbilt, Phys. Rev. B 41, 7892 (1990). 

[14] J. P. Perdew, J. A. Chevary, S. H. Vosko, K. A. Jackson, 
M. R. Pederson, D. J. Singh and C. Fiolhais, Phys. Rev. 
B 46, 6671 (1992). 

[15] S. Da,tta,,Electronic Transport in Mesoscopic Systems 
Cambridge (1995). 

[16] K. S. Thygesen, M. V. Bollinger and K. W. Jacobsen, 
Phys. Rev. B 67, 115404 (2003). 

[17] Y. Xue, S. Datta and M. A. Ratner Chem. Phys. 281, 
151 (2002) 

[18] M. Brandbyge, J. L. Mozos, P. Ordejon, J. Taylor and 
K. Stokbro Phys. Rev. B 65, 165401 (2002). 

[19] Y. Meir and N. S. Wingreen, Phys. Rev. Lett. 68, 2512 
(1992). 

[20] D. Djukic and J. M. van Ruitenbeek, private communi- 
cation 

[21] I. Souza, N. Marzari and D. Vanderbilt, Phys. Rev. B 

65, 035109 (2001). 

[22] K. S. Thygesen and K. W. Jacobsen, 



5 



arXivxond-mat/0411086 
[23] D. M. Newns, Phys. Rev. 178, 1123 (1969). 
[24] A. Thielmann, M. H. Hettler, J. Konig and G. Schon 

Phys. Rev. B 68, 115105 (2003). 
[25] T. K. Ng and P. A. Lee Phys. Rev. Lett. 61, 1768 (1988). 
[26] P. Hyldgaard, S. Hershfield, J. H. Davies and 



J. W. Wilkins Ann. Phys. 236, 1 (1994). 
[27] B. Hammer and J. K. N0rskov, Surface Science 343, 211 
(1995). 

[28] J. K. N0rskov, A. Houm0ller, P. K. Johansson and B. I. 
Lundqvist, Phys. Rev. Lett. 46, 257 (1981). 



